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Abstract. The expansion of an initially unmagnetized planar rarefaction wave has 
recently been shown to trigger a thermal anisotropy-driven Weibel instability (TAWI), 
which can generate magnetic fields from noise levels. It is examined here if the TAWI 
can also grow in a curved rarefaction wave. The expansion of an initially unmagnetized 
circular plasma cloud, which consists of protons and hot electrons, into a vacuum is 
modelled for this purpose with a two-dimensional particle-in-cell (PIC) simulation. It 
is shown that the momentum transfer from the electrons to the radially accelerating 
protons can indeed trigger a TAWI. Radial current channels form and the aperiodic 
growth of a magnetowave is observed, which has a magnetic field that is oriented 
orthogonal to the simulation plane. The induced electric field implies that the electron 
density gradient is no longer parallel to the electric field. Evidence is presented here for 
that this electric field modification triggers a second magnetic instability, which results 
in a rotational low-frequency magnetowave. The relevance of the TAWI is discussed 
for the growth of small-scale magnetic fields in astrophysical environments, which are 
needed to explain the electromagnetic emissions by astrophysical jets. It is outlined 
how this instability could be examined experimentally. 
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1. Introduction 

Energetic electromagnetic radiation is emitted by a wide range of astrophysical objects, 
such as supernova remnants (SNRs) [1], microquasars |2] and gamma ray bursts (GRBs) 
[3]. This radiation is at least partially attributed to the synchrotron- or synchrotron 
jitter emissions [1] of ultra-relativistic electrons, which gyrate in strong magnetic fields. 
Such extreme conditions do not exist in the interstellar medium and probably not in 
astrophysical jets or in SNR blast shells, if their plasma would be in a steady state. 
Efficient electron acceleration and magnetic field amplification mechanisms must be at 
work and they are driven primarily by collisionless plasma shocks [5]. 

Through their interaction with the electromagnetic fields in the shock transition 
layer, the electrons may gain enough energy to be injected into the diffusive shock 
acceleration process [6]. The latter is probably capable of accelerating charged particles 
to the highest observed cosmic ray energies and it should supply a population of highly 
relativistic electrons. The apparent magnetization of SNR shocks [7] is attributed to 
cosmic ray-driven instabilities, which can operate on large scales [8]. Filamentation 
instabilities of counterstreaming beams of charged particles are thought to be the origin 
of the micro-scale magnetic fields, in particular of those in GRB jets |9l HO], although 
the lifetimes of these fields might be too short [HI |12] . 

Filamentation or beam-Weibel instabilities develop in an initially spatially uniform 
system of two counter-streaming particle beams. Typically one examines this instability 
in the reference frame, in which the net charge and current vanish. The beams consist 
of electrons and ions and they may contain a sizeable fraction of positrons in GRB 
jets. The filamentation instability transforms the spatially uniform beams into a 
system of spatially separated current fiux tubes, which are enwrapped and separated by 
electromagnetic fields [13-17]. Observations of the filamentation instability in laboratory 
settings and the therefrom resulting electromagnetic radiation exist [18] and they are 
likely to occur close to shocks [13] in astrophysical settings. 

The Weibel instability [19] is another magnetic instability. It is driven by a thermal 
anisotropy of a single electron species. Weibel considered in his original work a spatially 
uniform bi-Maxwellian electron velocity distribution. The electrons were cooler along 
one direction than along the other two. The Weibel instability gives rise to the aperiodic 
growth of magnetowaves with a wavevector that is parallel to the cool direction [20] and 
with wavelengths that are comparable to an electron skin depth. The Weibel instability 
has been widely examined in the past both analytically and numerically [21-30]. Its 
importance for the generation of magnetic fields in jet outfiows and on a cosmological 
scale has been pointed out [31]. It has, however, previously been unclear how a strong 
thermal anisotropy could form and survive over spatial scales that exceed the wavelength 
of the modes, which are destabilized by the Weibel instability in its original form |19] . 

Such a process has been identified recently [32] . A plasma density gradient together 
with the high electron mobility gives rise to an ambipolar electrostatic field, which 
confines the electrons and accelerates the ions. The electric field vector is antiparallel 
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to the density gradient. The ion acceleration launches a rarefaction wave [33-38]. The 
electrons, which move from the dense plasma to the front of the rarefaction wave, are 
slowed down and reflected inelastically by the electric field of the expanding wave, 
through which they provide the energy for the ion acceleration. The modulus of the 
electron momentum along the density gradient is reduced by these processes, while both 
orthogonal momentum components are unchanged. The thermal anisotropy can give rise 
to the Weibel instability. In what follows we refer to it as the thermal anisotropy-driven 
Weibel instability (TAWI). Particle-in-cell (PIC) simulations [32] have shown that the 
TAWI can magnetize a planar plasma cloud that expands into a vacuum. A shock 
develops, if the rarefaction wave expands into a dilute ambient medium. The TAWI is 
suppressed if a shock is limiting the extent of the rarefaction wave to a length that is 
shorter than the wavelength of its modes [37l [38] . 

This process, which requires ion density gradients and hot electrons, can be 
important for energetic astrophysical plasma. The internal shock model for the prompt 
GRB emissions assumes that the jet consists of plasma clouds with varying densities 
and mildly relativistic relative speeds [3]. We can also assume that the ion density 
close to SNR shocks will not be spatially uniform. The bulk electrons in the GRB jets 
are thought to be relativistically hot and the bulk electrons close to SNR shocks have 
keV temperatures. The magnetic energy generated by the TAWI can reach a sizeable 
fraction of the electron thermal energy, which could result in strong magnetic fields. 

The TAWI has previously been examined numerically in a two-dimensional planar 
geometry. A plasma slab is centred in these studies in the simulation box and confined by 
two parallel boundaries, across which the density jumps from the constant high density 
within the slab to a low density or a vacuum outside of it. The plasma density is constant 
along the second simulation direction and the boundary conditions are periodic. The 
rarefaction waves launched at both boundaries move into opposite directions, causing 
an expansion of the slab plasma. The density gradients of both rarefaction waves are 
antiparallel. An electron can cross the slab many times and its momentum component 
along the density gradients remains decoupled from both orthogonal ones. Multiple 
electron reflections, which reduce only the electron momentum component along the 
plasma density gradient, result in an increasingly strong thermal anisotropy, which may 
exaggerate the strength of the TAWI. 

We examine here with a PIC simulation the expansion of an initially unmagnetized 
circular plasma cloud. The electron momentum can in this case not be subdivided into 
a component that is always parallel to the density gradient and into two components 
orthogonal to it, unless the electron moves exactly in the radial direction. This case is 
unlikely. Consecutive electron reflections by the electric field of the rarefaction wave will 
diminish different momentum components and the anisotropy between the temperature 
in the radial and azimuthal directions is not boosted by multiple reflections. The purpose 
of our simulation study is to assess if the slowdown of the electrons in the sheath potential 
alone is still strong enough to trigger the TAWI. 

Our results are as follows. The ambipolar electric field accelerates the protons in the 
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simulations to a few percent of c. A TAWI develops within the density gradient of the 
rarefaction wave. This instability is similar to that observed in Ref . |32] . It drives radial 
current channels, which yield the growth of a transverse magnetic (TM) field component. 
The growth of the TM wave is followed by that of a secondary TE magnetowave. The 
rotational polarization of this TE wave is concluded from the simultaneous presence 
of in-plane magnetic vortices and out-of-plane magnetic fields that interfere with those 
driven by the TAWI. The combined magnetic field of both waves breaks the circular 
symmetry of the rarefaction wave and scatters rapidly the electrons within, thereby 
reducing their thermal anisotropy. The instability saturates in our simulation, when the 
total magnetic energy density is about 1% of the electron thermal energy. 

The TE and TM magnetowaves develop in the same radial interval, which suggests 
their connection. We propose that the superposition of the electric field, which is induced 
by the magnetic field of the TAWI, with the radial electric field of the rarefaction wave 
triggers the growth of the TE wave by the misalignment of the plasma density gradient 
and the electric field vector. Such a system is unstable [391 HO] and it results in an 
instability similar to a thermoelectric instability. 

The structure of this paper is as follows. Section 2 describes the PIC simulation code 
and the initial conditions. Our results are presented in section 3. Section 4 summarizes 
our findings and it discusses how these could be verified experimentally. 



2. The simulation code and the initial conditions 

The 2D3V (resolves two spatial and three momentum dimensions) particle-in-cell 
simulation code we employ here is based on the virtual particle-mesh numerical scheme 
|41j . Faraday's and Ampere's laws 

VxE^-f (1) 

V X B = /ioJ + /ioeo^ (2) 

are solved on a grid and evolve the electromagnetic fields in time, while V ■ E = p/eo 
and V ■ B = are fulfilled as constraints. The code solves all three field components as 
a function of x = {x,y). The relativistic force equation 

^ = g,(E[x,]+v,xB[x,]) (3) 

is used to update the momentum of the i^^ computational particle (CP) of species j, 
with charge qj and mass rrij, that is located at the position Xj. The electric E and 
magnetic B fields are connected with the CPs as follows. The micro-current of each CP 
is interpolated to the grid. The sum of the interpolated micro-currents of all CPs gives 
the macroscopic current J, which is used to update the electromagnetic fields through 
Eq. [21 The new fields are interpolated back to the positions of the individual CPs to 
update their momentum in time through Eq. El The position of each CP is updated 
with diLi/dt = {vi^x,Vi^y). A detailed review of the PIC method is given elsewhere |i2] . 
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Figure 1. The 2D simulation plane covers the interval —L/2 < x,y < L/2. A circular 
plasma cloud has its center at a;,y = and the radius rw — L/5A. It is filled with 
electrons and protons with the density uq. The boundary conditions are periodic for 
the fields and for the particles. We will in some cases exploit the circular symmetry 
and use the circular coordinates (r, p) instead of (x, y). 



The plasma cloud with radius is surrounded by vacuum and it consists of 
electrons with the number density no and protons with the same density, which are 
initially distributed uniformly across the cloud. Figure [1] illustrates the setup. The 
electron's relativistic Maxwellian velocity distribution has a temperature of 32 keV or a 
thermal speed Ve = c/4. The proton thermal speed in both simulations is fp = 3.1 x 10^ 
ms~^ and corresponds to a thermal energy of 10 eV. The mean speeds of the electrons 
and of the protons vanish initially. The electrons and protons are approximated each 
by ~ 2500 CPs per cell. The electron mass is mg and the proton mass rup = 1836me. 

1 /2 

The Xe = c/oop is the initial electron skin depth and Up = (noe^/meeo) is the electron 
plasma frequency. The time is expressed in units of The time step is = 0.022. 
The simulation time is Tp = 880. The cloud radius = 14.2Ae. An electron with 
the thermal speed Ve will cross the cloud diameter 2rw about 8 times during t = Tp, 
resulting in multiple bounces at the electrostatic sheath field. The simulation box with 
periodic boundary conditions in all directions uses 1700 x 1700 cells for a square domain 
with the side length L = 75.5Ae. 

3. The simulation results 

We examine the plasma phase space distributions and the electromagnetic field 
distributions, the loss of electron thermal energy to the protons and the thermal 
anisotropy. We consider in more detail the simulation times Ti = 5.5, T2 = 550 and 
T3 = 750. This section concludes with a discussion of the time-evolution of the plasma 
quantities and effects due to the periodic boundary conditions. 
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Figure 2. (Color online) The field distributions at the time Ti. (a) shows the modulus 
of the in-plane electric field Ep. (b) displays the in-plane magnetic field Bp and (c) 
the orthogonal field \Bz\. A strong electric field points radially outwards and outlines 
the cloud's perimeter. The magnetic fields are at noise levels. Overplotted is a black 
circle with r = rw- 



We transform the distributions from a representation in (x, 2/)-space into one in 
(r, p)-space, where r = is the center of the cloud and of the simulation box. We 
exploit the approximate circular symmetry and average these quantities over the azimuth 
angle p. This average will become less accurate at late times, when B breaks this 
symmetry. The electron number density fie(r) = Jne{r,p)rdp. The phase space 
densities fe,p{r, p,Vr,Vp,Vz) are integrated as fe,p{r,Vr) = J f^^prdpdvpdvz, where 
e and p stand for electrons and protons. The phase space densities at all times are 
normalized to the maximum of the respective distribution at t = 0. The azimuthally 
averaged radial thermal energy Kr{r) = K^^ J^rdp J fev'^dv and azimuthal thermal 
energy Kp{r) = Kq^ J rdp J feVpdv define the thermal anisotropy A{r) = Kr{r) / Kp{r). 
The initial electron thermal energy per degree of freedom is Kq. The in-plane electric 
field Ep = \Ex + iEy\ and the orthogonal E^ field are expressed in units of UpCnie/e, 
while the magnetic field components Bp = {Br^ + iBy\ and B^ are expressed in units 
of UpiTie/e. The fields are azimuthally averaged as (F^) = (27rr)"^ F^rrfp, where F 
serves as a place holder for Ep,Ez, Bp and Bz- 

3.1. Early time Tl: Proton acceleration 

Their mobility implies that some electrons can escape from the cloud and build up a 
negatively charged sheath outside of the now positively charged cloud. An ambipolar 
electrostatic field develops that points along the surface normal. Its modulus is displayed 
in Fig. |2]for the time Ti. The electric field outlines the cloud radius and has a thickness 
comparable to Ae- The magnetic Bp and Bz components show noise that is strongest 
for r < rw and weaker in the rarefaction wave. The noise originates from statistical 
fluctuations of the current density. The initial plasma evolution is dominated by the 
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Figure 3. (Color online) The azimuthally averaged radial phase space density at the 
time Ti normalized to its initial value: Panels (a) and (b) show /e and /p, respectively. 
The color scale is 10-logarithmic and rw is the initial cloud radius. 



electrostatic fields because their amplitude is much larger than the amplitude of the 
magnetic fields and because the latter are incoherent. 

The protons are accelerated radially outwards by the electric field, which confines 
at the same time the electrons. This is confirmed by Fig. [31 The electrons lose speed 
in the radial direction and their density decreases drastically as they cross the strong 
in-plane electric field Ep aX r = r^. This field is not strong enough to hold back all 
electrons and some escape into the vacuum. The fastest electrons have expanded to 
a radius r ^ 1.3rw, while the proton expansion is still negligible. The protons with 
r ^ are accelerated by Ep and some reach already at the time Ti a peak speed 
that exceeds their initial thermal speed Vp by an order of magnitude. The protons with 
r < Tyj still show their initial velocity distribution. Their low thermal speed implies, 
that they could not yet interact with the surface electric field. 

3.2. Intermediate time T2: The TAWI 

Figure H] shows the particle distributions at the time T2. The protons have been 
accelerated by the ambipolar electric field to about 200fp or to 2% of the speed of 
light c. The protons form a narrow beam with a speed, which grows approximately 
linearly with r for r > O.Trvi/, and a density that decreases with an increasing radius. 
This high speed has been reached by the protons, because there is no shock that limits 
their acceleration. A shock would form, if the rarefaction wave would expand into an 
ambient plasma. The electron distribution has been altered significantly at this time. 
The peak speed of the electron's bulk population in the cloud's core, which has shrunk 
to r < 0.7riy at this time, has been decreased compared to that in Fig. EJ^a), which 
is indicative of an energy loss. This energy loss can be attributed to the electron's 
inelastic reflection by the ambipolar electric field of the expanding rarefaction wave. 
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Figure 4. (Color online) The azimuthally averaged radial phase space density at the 
time T2 normalized to its initial value: Panels (a) and (b) show /e and /p, respectively. 
The color scale is 10-logarithmic and rw is the initial cloud radius. 
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Figure 5. (Color online) The field distributions at I—T2: Panel (a) shows the modulus 
of the in- plane electric field Ep and (b) shows that of the in-plane magnetic field Bp. 
Panel (c) shows the modulus \Bz\. The black circles with radius rw outline the initial 
cloud perimeter. 



The electrons fill up the entire simulation box, although their phase space density close 
to the boundary is 4 orders of magnitude less than the peak density within r < 0.7rvi/. 

The erosion of the cloud implies, that the initial distribution of Ep in Fig. [2|^a) 
could not have been upheld. This is confirmed by Fig. O which shows the Bp and 
the modulus of Bz at the time T2. The electric field maximum is located at r 1.6riy, 
which is at the tip of the expanding protons in Fig. Hl^b). The electric field distribution 
expands to lower r. It thus continues to accelerate the protons and sustains the increase 
of the proton mean speed with the radius. The electric fields aXr rw are a consequence 
of the plasma density gradient in this interval and of the higher electron mobility, which 
implies a space charge. Figure [5]^a) reveals, that the electric field in the interval r ~ rw 
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Figure 6. Panel (a) shows the modulus of the in-plane electric field Ep{x,y). The 
contour lines correspond to the values -0.01, -0.005, 0.005 and 0.01 of the field. 
The Bz field has maxima in the upper right and lower left quadrant and minima in the 
other two quadrants. Panel (b) compares the same contour lines of Bz with the contour 
lines 2.7, 2.9, 3.1 and 3.3 of the 10-logarithmic electron number density. The density 
increases with decreasing radius. The number density is expressed in computational 
particles per cell and logj^g ng = 3.4. 



does not show an exact circular symmetry. The amphtude is a function of p. 

This break of the circular symmetry and its origin become evident in Fig. |6l which 
shows the amplitude modulus of Ep and the electron density distribution in the relevant 
spatial interval and compares both quantities to B^. The amplitude of the azimuthal 
oscillation of Ep at r ^ rw is about 10% of the mean value. Figure EJ^a) suggests a 
correlation of Ep with Bz for two reasons. Firstly the azimuthal oscillation of Ep and 
the strong Bz are observed at approximately the same radius. Secondly the Ep field 
shows two amplitude minima and two maxima, which are shifted by 45° relative to 
those of Bz- We find one maximum of Ep aX x ^ rw and ?/ ~ (p ~ 0°), while one 
maximum of Bz is located at x,y ^ rw/V^ (p ~ 45°). 

The circular symmetry of Ep can, in principle, be broken by an electrostatic drift 
wave along the azimuthal direction. They could be generated by the lower-hybrid 
drift instability |l3l Hi], which is caused by gradients in the plasma density and in 
the magnetic field amplitude. A second possibility is the Simon-Hoh instability that 
develops, if a radial electrostatic field and a spatially uniform orthogonal magnetic field 
are present |l5l US] . The electron E x B drift destabilizes the plasma. The electric field 
driving the Simon-Hoh instability would be the Ep and the orthogonal magnetic field 
would be Bz- However, the characteristic wave lengths of the waves driven by both 
instabilities would be much less than that of the modulation of Ep observed in Fig. 
[6l The lower hybrid instability destabilizes wavelengths, which are shorter than the 
characteristic length scale of the magnetic field gradient and the Simon-Hoh instability 
requires that Bz be uniform over a scale much larger than the wave length of the unstable 
wave. Figure [6](b) furthermore demonstrates that the electron density distribution is 
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Figure 7. The azimuthally averaged plasma configuration at t—T2: The (dashed) 
curve rig corresponds to the radially averaged electron number density fii,/r, which 
is normalized to riQ. The curve ujp equals {fie /nor) ' and denotes the local plasma 
frequency. The curve Eq is the total magnetic energy density r{Bp) + r{Bl). The 
thermal anisotropy A = Kr{r)/ Kp(r) is the second dashed curve. 

still circularly symmetric. The azimuthal oscillation of Ep is thus not correlated with a 
charge density modulation and it is probably not an electrostatic wave. 

We can understand the modulation of Ep, if we take into account the in-plane 
current ( Jr, Jy) and electric field {E^, Ey) that sustain the localized B^. We neglect the 
weak Bp in Fig. |5]^b). Ampere's law (Eq. |2]) connects B^ with the in-plane electric field 
and the in-plane current. Consider the right part of the simulation box x > in Fig. 
El^a). The contour lines with ?/ > denote a maximum of B^ and those with ?/ < a 
minimum. At x ~ and y we find with the right hand rule that the electric field 
induced by V x (0, 0, Bz) gives dtE^ > and dfE^ ^ \dtEy\. The electric field amplitude 
should grow to a sizeable amplitude before a balance between the in-plane current and 
Bz cancels any further growth of Ep. The superposition of the ambipolar electric field 
from the rarefaction wave and the practically parallel induced iJ^^ at x ~ rw and ?/ ~ 
imply that Ep is strong. The induced electric field will also amplify the ambipolar 
electric field at x ~ —rw and y ~ 0, because both fields reverse their sign. The induced 
electric field weakens the ambipolar electric field at x ~ and |?/| ~ r^. A magnetically 
induced electric field is not tied to charge density modulations, which would explain 
why the circular symmetry is maintained for the electron density in Fig. EJ^b). 

We turn now to the identification of the source mechanism of the magnetic B^ 
component in Fig. [5]^c). Figure [7] compares the azimuthally averaged energy density 
r{{Bp) + {BD) of the magnetic field, the thermal anisotropy A{r) = Kr{r)/Kp{r) and 
the azimuthally averaged electron number density fie/r. The division of fie by r ensures 
that a uniform density within the cloud results in a constant value of Ue/r and the 
multiplication of the azimuthally averaged energy density by r gives us a curve that is 
proportional to the magnetic energy. These quantities are computed at the time T2. 
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Figure 8. (Color online) The azimuthally averaged radial phase space density at the 
time T3 normalized to its initial value: Panels (a) and (b) show /e and /p, respectively. 
The color scale is 10-logarithmic and rw is the initial cloud radius. 

The onset of the rarefaction wave is found at r ^ Q.Qr^ in Fig. [TJ where the electron 
density and plasma frequency start to decrease. A sidelobe of the magnetic energy Eb 
is found within this radius, but most energy is concentrated within 0.7 < r/rw < 1.1. 
The thermal anisotropy curve is a steadily decreasing function up to r lAryy and it is 
below 1 everywhere. The mean kinetic energy of the electrons in the radial direction is 
thus lower than their azimuthal one, which we attribute to the radial proton acceleration. 
The TAWI destabilizes waves with wave vectors along the cool direction. Figure ^c) 
shows that Bz oscillates fastest in the radial direction. One oscillation is visible with 
the wavelength ^ O.Sr^y. The magnetic component performs two oscillations in the 
azimuthal direction, each with a wavelength ^ nrw- A long oscillation orthogonal to 
the density gradient has also been observed in simulations of a planar rarefaction wave 
[32], which is probably a consequence of the density gradient that has not been taken 
into account by the existing analytic studies. 

The simulation data suggests that the large anisotropy observed in Fig. [7] triggers 
the TAWI. The Eb peaks at r ^ 0.85rw, where A = 0.7. One may think that Eb 
should peak in the interval, where A is lowest. However, both the growth rate of the 
TAWI in units of s~^ and the magnetic field amplitude the plasma can support increase 
with Up. The latter can be concluded from the fact, that we can normalize the Vlasov- 
Maxwell set of equations with the help of Up. A higher Up results in a higher B for the 
same plasma configuration. The magnetic amplitude decreases to noise levels at about 
r/rw ~ It is thus confined to within the rarefaction wave and it does not extend 
into the surrounding electron plasma at this time. 

3.3. Final time T3: Secondary instability 

The electrons have spread out in significant numbers up to r ~ 2rw at the time T3. The 
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Figure 9. (Color online) The field distributions at t=r3: Panel (a) shows the in-plane 
electric field Ep and (b) shows the in-plane magnetic field Bp. Panel (c) shows the 
modulus \Bz\. The black circles with radius rw outline the initial cloud perimeter. 



radial electric field has accelerated the protons to about 260fp or 2.6% of c at the time 
T3, as we can see from Fig. |8]^b). The protons form a cool beam with one well defined 
mean speed for any r < l.Srvi/- This compact distribution justifies its approximation 
by a single cold fluid [31]. The mean speed of the protons increases linearly within 
0.7 < r/rw < 1.6. 

Figure |9] shows -Ep, Bp and \Bz\ at the time T3. The in-plane electric field has 
weakened while the magnetic fields have grown further. The plasma density gradient 
and the ambipolar electric field in a rarefaction wave decrease with time. The infiuence 
on the particle dynamics of the magnetic field and its induced electric field thus become 
more important. The circular symmetry of Ep imposed by the rarefaction wave is 
no longer dominant in Fig. [9]^a). The amplitude of the in-plane Bp in Fig. [9]^b) 
has increased by an order of magnitude compared to that in Fig. [5](b) and it is now 
comparable to the amplitude of B^. The distribution of in Fig. Ws) does no longer 
show four almost equally important peaks as it did in Fig. [5t^c). 

We would expect such a pattern, if the wave responsible for the Bp is driven by 
a separate secondary instability and if it results in a circularly or elliptically polarized 
magnetowave. An instability is discussed in the next subsection that can generate such 
a wave. The Bp field oscillates much faster in the radial than in the azimuthal direction 
in Fig. [U^b) and its wave vector is thus practically aligned with the radius vector. We 
call this magnetowave wave 2. An elliptically polarized wave will have a B^ component 
that is phase shifted along r with respect to its Bp component. The B^ component of 
wave 2 would interfere with the B^ component of the wave that is driven by the TAWI 
(wave 1). Wave 2 performs one oscillation in the azimuthal direction according to Fig. 
|9] and we may assume that the same is true for its B^ component. The B^ field of 
wave 1 oscillates twice along p. So we expect that wave 2 interferes constructively with 
two half-oscillations of wave 1 and that it interferes destructively with the other two. 
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Interference between both waves could thus explain the pattern observed in Fig. [Hl^b). 
3.4- Time evolution 

Figure [7] has demonstrated that the magnetic field energy density is increased in 
the radial interval, in which we find a strong thermal anisotropy in the electron 
distribution. The anisotropy is caused by the loss of electron thermal energy to the 
radially accelerating protons. The supplementary movie 1 animates in time this energy 
exchange. It shows the electron phase space distribution in (a) and the proton phase 
space distribution in (b), which are displayed in the form of snapshots and discussed in 
Figs. 131 m and |8l Movie 1 confirms that electrons, which cross the boundaries on one 
side, re-enter on the opposite side, as we expect from the periodic boundary conditions 
of the simulation. The fastest electrons that escaped from the plasma cloud reach the 
boundaries at t ^ 50 and they encounter the electrons that crossed the boundary from 
the other direction. Note that movie 1 (a) shows the velocity modulus and does thus not 
reveal that both electron populations move into opposite directions. The electron two- 
stream configuration can result in instabilities. However, the two-stream configuration 
is rapidly destroyed by the arrival of the slower electrons at the boundaries. 

The supplementary movie 2 demonstrates that the life-time of the two-stream 
configuration is too short and the electron density too low to drive detectable instabilities 
at t ~ 50. It animates in time the 10-logarithmic electron density in the simulation box, 
which is averaged over blocks of 4 x 4 cells and given in particles per cell. An average 
density of less than one computational particle per cell is found close to the boundaries. 
The local plasma frequency is ~ Wp/lOO, which implies a low growth rate for any beam 
instability [161 [U] • Movie 2 reveals faint electron structures at a time t 200 that fiow 
parallel to the x-axis and y-axis. They are probably caused by an instability between 
the electrons of the plasma cloud and the electrons, which return from the boundary 
and are accelerated by the ambipolar electric field when they re-enter the cloud 

Figure Uni quantifies both the energy loss of the electrons to the accelerating protons 
and it follows the thermal anisotropy in time. Figure [TOT a) is a measure for the average 
kinetic energy per electron, which corresponds to their velocity in the simulation plane. 
This energy is obtained by summing Kr{r) and Kp{r) and by dividing the result by 
2KQhe/nQ. Kq is the initial thermal energy contributed by each velocity component and 
the division by fie/ no gives us a measure of the kinetic energy per electron. The thermal 
anisotropy A{r,t) is displayed in Fig. [TOT b). 

The electrons within the cloud lose about 30% of their kinetic energy to the 
accelerating protons. The electron kinetic energy is apparently increased outside of 
the rarefaction wave. This is a selection process, which can be seen in the movie 1. 
The bulk of the slow electrons is contained by the rarefaction wave, while the most 
energetic electrons, which have a kinetic energy that far exceeds the potential energy of 
the rarefaction wave, can escape. 

The thermal anisotropy is negligible for the radial interval without a plasma density 
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Figure 10. (Color online) Panel (a) displays the time evolution of the electron 
thermal energy nf){Kr{r) + Kp{r))/2Kohe and panel (b) the thermal anisotropy 
A{r) = Kr{r)/ Kpir). The colour scale is linear. The horizontal black line corresponds 
to the time T2 



gradient, e.g. for r < CSr^i/ and t = T2. A{r) reveals that the electrons in the interval 
with the proton density gradient lose more energy in the radial than in the azimuthal 
direction. The radial electron energy is less than half of the azimuthal energy over an 
extended radial interval. The interval with a value of A{r) < 0.6 is largest at the time 
T2 in Fig. [TOT b). when it covers 0.9 < r/rw < 1.3. The saturates at around this 
time. The TAWI starts to grow at t ^ 200, when the anisotropy is sufficiently strong. 

The saturation is demonstrated by Fig. [TTl which compares the time-evolution 
of the mean amphtude of the azimuthally averaged fields (Ep)^'^, {E'^)^'^, (-8^)°'^ and 
{Bl)^'^ . The supplementary movie 3 animates in time the spatial distributions of the 
in-plane electric field Ep (a), the in-plane magnetic field Bp (b) and the Bz- The in-plane 
electric field Ep in Fig. [TlT a) peaks initially at r = rw- The build-up of this electrostatic 
field launches electrostatic waves into the cloud, which yields the amplitude peaks for 
r < rw- The electric field that expands to larger r > and reaches r ^ 2rw at 
t = 50 is induced by the electrons that escape from the cloud, e.g. those with r > ryy 
in Fig. [3t^a). The electron current is not balanced by an ionic current and a radial 
electric field grows. The density of the escaped electrons decreases as they spread over 
the simulation box and the electric amplitude eventually drops below the threshold of 
Fig. [TlT a). The position associated with the maximum of Ep accelerates in time. This 
maximum coincides with the tip of the rarefaction wave. The electric field with the 
uniform amplitude of 10~^ behind the tip of the ion front roughly outlines the extent of 
the rarefaction wave. It expands to lower r < as the cloud is progressively eroded. 
The electric field component E^ in Fig. ITlT b) remains at noise levels, except for a 
peak at if: ~ T2. This peak coincides spatially and temporally with the growth of the 
rotational magnetowave, which is responsible for the -Bp 7^ in Fig. ITTT c). It is related 
to the displacement current dE^/dt connected to Bp. 
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Figure 11. (Color online) The azimuthally averaged electromagnetic fields (-E^) 
(a), {ED^-"^ (b), {Bl)°-^ (c) and (S^)"'^ (d). The color scale is lO-logarithmic. The 
black horizontal line corresponds to the time T2. 

A comparison of the Figs. [11] (c) and (d) confirms the observation from the 
supplementary movie 3 that starts to grow to a large amplitude well before Bp. 
The latter grows fastest at around the time t ~ T2, when B^ saturates. The presence of 
a strong B^ is apparently necessary for the growth of i?p, which supports our previous 
claim that Bp is driven by a secondary instability. Figure ITT] reveals that Bz is confined 
to within the rarefaction wave, while the magnetic field distribution Bp reaches out 
into the dilute electron plasma outside of the rarefaction wave. It eventually reaches 
the boundaries and effects due to the periodic boundary conditions can no longer be 
neglected at the final time Tp. The development of the strong magnetic fields aX r ry/ 
coincide with a balancing of and Kp in Fig. [TUT b). This change in A is partially due 
to the electron deflection by B, which will contribute to an equilibration of and Kp. 

We want to determine the source mechanism of the in-plane magnetic field Bp. 
The ambipolar electric field due to the rarefaction wave and the electric field induced 
by Bz are both confined to the simulation plane. The electric fields and B^ can excert 
a force on the electrons only in this plane and no current orthogonal to the simulation 
plane can be generated. Since we have no initial beams moving along the z-axis, we 
can rule out magnetic beam instabilities [lEl [171 121]. One possibility to generate an 
electron flow along z is a coupling of the electron velocity components through the 
Lorentz factor. An electron acceleration in the simulation plane would change Vz- The 
fastest electrons in Fig. Il](a) have the mildly relativistic speed v ~ 0.7c that gives a 
Lorentz factor F ^ 1.4. However, their number density is low and we could not expect 
them to drive strong in-plane magnetic fields. A possible mechanism for the growth of 
Bp is the electric field induced by Bz- The induced field is strong when B^ saturates 
and it changes the direction of the in-plane electric field vector by up to 5% relative to 
the radius vector. The plasma density is according to Fig. |6]J^b) circularly symmetric 
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at t = T2. A misalignment of the in-plane electric field vector and the density gradient 
can give rise to elliptically polarized magnetowaves [391 

4. Discussion 

We have examined here the expansion of an initially unmagnetized circular plasma 
cloud into a vacuum by means of a 2D3V PIC simulation. The purpose of our study 
has been to test, if a thermal anisotropy-driven Weibel instability (TAWI) can also 
develop in curved rarefaction waves. Previous studies considered planar plasma slabs 
that expanded under their own thermal pressure [32] and observed the growth of strong 
magnetic fields in the rarefaction wave, unless a shock suppressed the TAWI [371 138] . 

A thermally anisotropic electron distribution develops due to two reasons. Firstly 
the slow-down of the electrons in the ambipolar electric field implies that their thermal 
pressure is decreased along the density gradient. This mechanism occurs in planar 
and curved rarefaction waves. Secondly, the inelastic reflection of an electron by 
the expanding rarefaction wave implies that it looses momentum along the electric 
field direction. This momentum is transferred to the accelerating ions. The initial 
plasma distribution used in the Refs. [321 |371 [3H] is a plasma slab of infinite extent 
in one direction, which has a finite width in the second direction. This plasma slab is 
immersed in a vacuum or in a low-density plasma. Both plasma boundaries are parallel. 
This implies that every reflection of the electron will decrease the same momentum 
component, namely the component that is aligned with the antiparallel density gradients 
of both rarefaction waves. The anisotropy increases with every electron reflection. This 
mechanism can not occur in a circular rarefaction wave, since here consecutive electron 
reflections tend to decrease different momentum components. Our study can thus test 
if consecutive electron reflections are necessary to achieve a TAWI. 

If the slowdown of the electrons in the ambipolar electric field of a plasma with a 
spatially non-uniform density is enough to trigger the TAWI, then this instability may 
play a role in the magnetization of astrophysical plasma [31] besides the filamentation 
instability [9] and cosmic ray-driven instabilities [8]. The magnetic energy that can be 
generated by the TAWI is a fraction of the electron thermal energy, here about 1%. This 
energy density is well below that of the filamentation instability, which can reach about 
10-20% of the kinetic energy density of leptonic beams if they are cold. However, the 
TAWI requires only hot electrons and a density gradient to grow, while the filamentation 
instability is triggered by relativistic particle beams that we may find only close to 
plasma shocks. The TAWI could thus magnetize much larger intervals in energetic 
astrophysical jets, such as GRB jets, which are thought to have a nonuniform plasma 
density and relativistically hot electrons. The TAWI may also generate a seed magnetic 
field in the turbulent plasma close to SNR shocks and speed up the development of the 
much more powerful cosmic ray-driven instabilities. 

Our findings follows. The electron expansion results in a strong radial 

electrostatic field. It counteracts a further electron expansion and accelerates the protons 
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into the radial direction. A rarefaction wave emerges, which is characterised by an 
increase of the proton mean speed and by a decrease of the proton density with an 
increasing radius. The expansion of a rarefaction wave is self-similar [331 EH |35] . The 
absence of an ambient medium has prevented the formation of a shock and the protons 
could reach a speed exceeding 250 times their thermal speed. 

The fields within the rarefaction wave have initially been electrostatic. Electrons 
lose momentum to the accelerating protons only in the radial direction, which implies 
that their radial thermal pressure decreases more quickly than the azimuthal one. A 
thermal anisotropy-driven Weibel instability (TAWI) [19j grows once the radial extent 
of the rarefaction wave is sufficient to accommodate a wavelength of the unstable modes 
[32] and once the thermal anisotropy between the radial and the azimuthal direction 
became sufficiently large. The TAWI results in the aperiodic growth of a transverse 
magnetic (TM) wave. This TM wave induces an electric field that is added to that of 
the electrostatic ambipolar electric field of the rarefaction wave. The total electric field 
is thus no longer parallel to the plasma density gradient. 

A misalignment between the electric field and the plasma density gradient can drive 
circularly or elliptically polarized magnetowaves [391 HO] • We propose that these waves 
are responsible for the in-plane magnetic field we observe in our simulation, which can 
not be explained in terms of the TAWI. The magnetic fields of this wave reach an 
amplitude that is comparable to the waves due to the TAWI. 

Our findings can be verified in an experiment, where the absorption of a laser pulse 
heats up the electrons in a thin wire [IE]. Some electrons will leave the now positively 
charged wire and build up an ambipolar electric field. This ambipolar field accelerates 
the ions through the transverse normal sheath acceleration (TNSA) mechanism [49l [50] 
and it results in the formation of a rarefaction wave. The ablation of a wire generates 
a cylindrical rarefaction wave. Our circular plasma cloud corresponds to an idealized 
cross section of the wire, as long as it is sufficiently far away from the laser impact point 
[18] . The cross section is idealized, because we do not consider ionization processes 
that require more elaborate numerical codes [511 ES] . The plasma dynamics outside of 
the wire should be comparable in the simulation and experiment. The strong magnetic 
fields that develop in the rarefaction wave should be detectable. It is thus possible 
to investigate systematically the TAWI in a controlled laboratory experiment and to 
determine the strength and life-time of the magnetic field. The interplay of the TAWI 
with the secondary magnetowave we have observed here and other instabilities, such 
as the Ray leigh- Taylor instability [531 El]; can also be examined. The experiment 
can reveal the topology of the magnetic fields in a three-dimensional setting, which 
is currently computationally too expensive. 
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